Physical regularization for the spin-1/2 Aharonov-Bohm problem in conical space 
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We examine the bound state and scattering problem of a spin-one-half particle undergone to an 
Aharonov-Bohm potential in a conical space in the nonrelativistic hmit. The crucial problem of 
the 5-function singularity coming from the Zeeman spin interaction with the magnetic flux tube is 
solved through the self-adjoint extension method. Using two different approaches already known in 
the literature, both based on the self-adjoint extension method, we obtain the self-adjoint extension 
parameter to the bound state and scattering scenarios in terms of the physics of the problem. It is 
shown that such a parameter is the same for both situations. The method is general and is suitable 
for any quantum system with a singular Hamihonian that has bound and scattering states. 
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Singularities are very common in quantum mechanics 
and already have a long history [1]. The first work with 
(5-like singularities was in the Kronig- Penny model [2] for 
the description of the band energy in solid-state physics. 
In addition, point interactions [3-5] have been of great 
interest in various branches of physics for their relevance 
as solvable models [6]. In the Aharonov-Bohm (AB) ef- 
fect [7] of spin-1/2 particles [8-10] a two-dimensional 6 
function appears as the mathematical description of the 
Zeeman interaction between the spin and the magnetic 
fiux tube [11, 12]. Hagen [9] argued that a ^- function con- 
tribution to the potential can not be neglected when the 
system has spin, having shown that changes in the am- 
plitude and scattering cross section arise when the spin 
of the particle is considered. Point interactions usually 
appear in quantum systems in the presence of topological 
defects. A simple but nontrivial example is the case of a 
cone rising from an effective geometry immersed in sev- 
eral physical systems, such as cosmic strings [13], defects 
in elastic media [14], defects in liquid crystals [15], and 
so on. In such systems, although the particle does not 
have access to the core (defect) region, its wave function 
and energy spectrum are truly influenced by it. 

Recently, a device was proposed that would detect 
microstresses in graphene [16] based on a scanning- 
tunneling-microscopy setup able to measure AB interfer- 
ences at the nanometer scale. In this setup a (5-function 
scattering potential was considered in the continuum 
limit [17]. In Ref. [18] it was considered a topological 
insulator nanowire with a magnetic field applied along 
its length, focusing on the AB conductance oscillations 
arising from the surface states. The Dirac Hamiltonian 
of this model takes into account the spinorial connection 
that allows us to incorporate topological defects (aris- 
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ing from a nontrivial conical geometry) through the met- 
ric. From these studies, such materials could be analyzed 
through theoretical models allowing to include point in- 
teractions able to reproduce AB-like effects. 

In quantum mechanics, singularities and pathological 
potentials, in general, are dealt with some kind of reg- 
ularization procedure. A common approach to ensure 
that the wave function in the presence of a singularity is 
square-integrable (and therefore might be associated to 
a bound state) is to force it to vanish on the singularity. 
More appropriately, an analysis based on the self-adjoint 
extension method [19], broadens the boundary condition 
possibilities that still give bound states. The physics of 
the problem determines which of these possibilities is the 
right one, leaving no ambiguities [8, 20]. This method 
has been applied by many authors, in particular, for AB- 
like systems [8, 12, 21-24]. However, the results obtained 
in these works present the most important results (e.g., 
energy spectrum, phase shift, S matrix) in terms of an 
arbitrary real parameter, the so called self-adjoint exten- 
sion parameter. 

In this article, we describe a general regularization 
procedure to obtain the self-adjoint extension parame- 
ter, based on the physics of the spin-1/2 AB system in 
(1 -|- 2) -dimensional conical space for both bound and 
scattering scenarios. We take as a starting point the 
works of Kay-Studer (KS) [25] and BuUa-Gesztezy (BG) 
[26], both based on the self-adjoint extension method. 

The topological defect considered here is a linear quan- 
tity that appears embedded in the metric system ds^ = 
dr"^ + a^r'^dip^ ^ where r > 0, < (/? < 27r, and a is the 
parameter which effectively introduces an angular excess 
or deficit, identified by 27r(l — a). The above metric has 
a conelike singularity at ?' = 0. In other words, the cur- 
vature tensor of this metric, considered as a distribution, 
is given by = R\ = R2 = 2'n{a — l)5{r)/a, where 
5{r) is the two-dimensional 5 function in fiat space [27]. 
This implies a two-dimensional conical singularity sym- 
metrical in the z axis, which characterizes it as a linear 
defect. 
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In order to study the dynamics of the particle in 
a nonflat spacetime, we should include the spin con- 
nection in the differential operator and define the re- 
spective Dirac matrices in this manifold. This system 
is governed by the modified Dirac equation in curved 
space [i^^{d^ + T^) - (77^7!^ - M]^{x) = 0, where q 
is the charge, M is mass of the particle, is a four- 

component spinorial wave function and is the spin 
connection. The only nonvanishing spin connection in 
this case is T^p ~ i{l — a)az/2, while the Dirac matrices 
are conveniently defined as a' = 7^7*, /? = 7" [28, 29]. 

The magnetic fiux tube in the background space de- 
scribed by the metric above considered is related [29] to 

the magnetic field sB = s(V x A) = ^^^z (where 
<j) ~ (j) / 2ti \s the flux parameter), while the vector poten- 
tial in the Coulomb gauge is A.^p = ^(p, with s ~ ±1 
being twice the spin projection parameter. The parame- 
ter s implies that the Dirac equation describes the planar 
motion (in the absence of the z coordinate) of the particle 
having only one projection of three-dimensional spin vec- 
tor. To examine the physical implications of these equa- 
tions, we consider their nonrelativistic limit. In this con- 



text, writing ?/' = {^^X) e 
equation is = idt^, with 



iMt 



the Schrodinger-Pauli 



H 



1 

2M 



-V„ - — 

I ar 



1-a 

c 

2ar 



qscj) S{r) 
2Ma~r~ 



(1) 



^ ^ -^-i^ is Laplacian op- 



where 

erator in the conical space, and (Ji — (dr, c,^, Cz) are the 
Pauli matrices in cylindrical coordinates. 

For this system the total angular momentum operator, 
J = — iVip -|- (72/2, commutes with the effective Hamilto- 
nian. So, the solution to the Schr"dinger-Pauli equation 
can be written in the form 
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with m = n + 1/2, n G Z. At the same time, the radial 
equation for fi(r) becomes 



where 



T~L ~ Ho + Z^short , 



^°-~2M 



1 d 
r dr 
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(5) 
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with j = ■^{ni — ^ ^ q<j> + ^^Y^)- The Hamiltonian in Eq. 
(4) governs the quantum dynamics of a spin-1/2 charged 
particle in the conical spacetime, with a magnetic held 



B along the z-axis, i.e., a spin-1/2 AB problem in the 
conical space. Let us consider a conical defect with a 
nucleus with radius ro , so it is suitable to write Z^short (r) 
as [9, 11] 



short 



(r) 



qs<j> S{r — ro) 



2Ma 
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and, at the end, the limit tq is taken. Although 
the functional structure of Wshort and Z^shoit arc quite 
different, as discussed in [9], we are free to use any form 
of potential provided that only the contribution of the 
form (6) is excluded. 

The operator 'Ho, with domain 2?('Ho)i is self-adjoint 
if 'D{'H\) = 'D{'Ho) and 'H\ = "Ho- For smooth functions, 
g e C(f'(R2) with 5(0) = 0, we should have Ug = Hog, 
and hence it is reasonable to interpret the Hamiltonian 
(4) as a self-adjoint extension of 'Ho\cg°(R'^\{o}) [30-32]. 
In order to proceed to the sclf-adjoint extensions of (5), 
wc decompose the Hilbert space = L^(R^) with re- 
spect to the angular momentum Jif = M'r ® , where 
= L2(M+,rdr) and = L'^{S^,dLp), with 5^ de- 

noting the unit sphere in K . The operator 



IS es- 



sentially self-adjoint in L'^iS^ ^ dip) [19] and we obtain the 
operator 'Ko in each angular momentum sector. Now, us- 
ing the unitary operator V : L'^{R~^,rdr) — > L^(IR+,d?-), 
given by iVg){r) = r^^^g{r), the operator Hq becomes 



ho = VHoV-^ = - 



r ^2 



2M 



dr^ 4 



1\ 1 
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which is essentially self-adjoint for |j| > 1, while for 
\j\ < 1 it admits a one-parameter family of self-adjoint 
extensions [19], Ho.Xj , where Xj is the self-adjoint exten- 
sion parameter. To characterize this family, we will use 
the KS [25] and the BG [26] approaches, both based in 
boundary conditions. 

In the KS approach, the boundary condition is a match 
of the logarithmic derivatives of the zero-energy solutions 
for Eq. (3) and the solutions for the problem Ho plus 
self-adjoint extension. In the BG approach, the bound- 
ary condition is a mathematical limit allowing divergent 
solutions of the Hamiltonian (5) at isolated points, pro- 
vided they remain square integrable. 

Now, the goal is to find the bound states for the Hamil- 
tonian (4). Following [25], we temporarily forget the S- 
function potential and find the boundary conditions al- 
lowed for Ho- But the sclf-adjoint extension provides 
infinity possible boundary conditions, so that it cannot 
give us the true physics of the problem. Nevertheless, 
once the physics at r = is known [8, 33, 34], it is pos- 
sible to determine any arbitrary parameter coming from 
the self-adjoint extension, so that it is possible to obtain 
a complete description of the problem. Since we have a 
singular point, we must guarantee that the Hamiltonian 
is self-adjoint in the region of motion. Note that even if 
Til ~ their domains could be different. 

We must find the deficiency subspaces, Af± , with di- 
mensions n+ and n_, respectively, which are called dc- 
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ficicncy indices of Hq [19]. A necessary and sufficient 
condition for Hq being essentially self-adjoint is that 
n_|- = n_ = 0. On the other hand, if = n_ > 1, 
then "Ho has an infinite number of self-adjoint extensions 
parametrized by unitary operators U : A/V — > Af- . 

Next, we substitute the problem in Eq. (3) by Hofg = 
£fg, with fg labeled by a parameter g which is related to 
the behavior of the wave function in the limit r — > tq. But 
we cannot impose any boundary condition (e.g. f ~ 
at r = 0) without discovering which boundary conditions 
are allowed to Hp- Then, from Eq. (5) we achieve the 
modified Bessel equation (k^ = —2M£, £ < 0) 



1 d 

r dr 



Mr) = 0. 



(9) 



Now, in order to find the full domain of Tin 



L^(IR+,rdr), we have to find its deficiency subspace. To 
do this, we solve the eigenvalue equation 



where Ho is given by Eq. (5). The only square-integrable 
functions that are solutions of Eq. (10) are the modified 
Bessel functions Ky^ {r^/^e), with e = 2iM. These func- 
tions are square integrable only in the range |j| < 1, for 
which Ho is not self- adjoint. The dimension of such de- 
ficiency subspace is (n_|-,77,_) = (1,1). Thus, 2?('Ho) in 
L'^{M.'^ ,rdr) is given by the set of functions [19] 

f,ir) = fiAr) + C [K\,\{rV~e) + e^^Ky\{rV~s)] , (11) 

where /i,j(r), with /ij(ro) = fi,j{ro) = (/ = df/dr), is 
the regular wave function when we do not have iYshort(?')- 
The last term in Eq. (11) gives the correct behavior for 
the wave function when r — ro- The parameter £i(mod27r) 
represents a choice for the boundary condition. As we 
shall see below, the physics of the problem determines 
such a parameter without ambiguity. In fact, g describes 
the coupling between Z^siiort('') and the wave function. 
Thus, it must be expressed in terms of a, the defect 
core radius rp and the effective angular momentum j. 
The next step is to find a fitting for g compatible with 
Wshort('')- In this sense, we write Eq. (3) for £ — 0, 
implying the zero-energy solution, Hfo = 0. Now, we 
require the continuity for the logarithmic derivative 



fo 



k 



(12) 



where fg{r) comes from Eq. (11). However, since tq w 0, 
the right-hand side of the Hamiltonian (12) is calculated 
using the asymptotic representation for Eq. (11) in the 
limit r — > 0. The left-hand side of Eq. (12) is achieved 
integrating the equation Hfo = 0, from to tq, which 
yields the parameter g in terms of the physics of the 
problem, i.e., the correct behavior of the wave functions 
for r — > rg. By solving Eq. (12) for £, we find the energy 



spectrum 




i/lil 



(13) 



Notice that there is no arbitrary parameters in the above 
equation. 

The above approach has the advantage of yielding the 
self-adjoint extension parameter in terms of the physics 
of the problem, but it is not appropriate for dealing 
with scattering problems. On the other hand, the BG 
method [26] is suitable to address both bound and scat- 
tering scenarios, with the disadvantage of allowing ar- 
bitrary self-adjoint extension parameters. Now, we ap- 
ply the BG approach to solve bound and scattering 
problems. By comparing the results of these two ap- 
proaches for bound states, the self-adjoint extension pa- 
rameter can be determined in terms of the physics of the 
problem. Here, all self-adjoint extensions of Ho,\ are 
parametrized by the boundary condition at the origin 



(10) [26] (go(r) =lim,^o+'^l^l <7(0) 



goir) = Xj 



lim 

r-yO+ rlJ 



gir) - go{r')- 



(14) 



The solutions for Hofij = k^fij 
can be written as (p = 2ikr) 



(fc2 = 2M£) for r 7^ 0, 



fiAr) 



Ml 



|j|,l + 2|j|,p) 



-Bje- 



p-^'^ iF,{l-\j\,l-2\j\,p), (15) 

where iFi(a,b, z) represents the confluent hypergconiet- 
ric function, and Aj , Bj are the coefficients of the regular 
and irregular solutions, respectively. By implementing 
Eq. (15) into the boundary condition (14), we derive the 
following relation between the coefficients Aj and Bj: 



(2z/c)-2blBj 1 



4(1 



— hm r 

IjI) r-^0+ 



2-2|j| 



(16) 



In the above equation, the coefficient of Bj diverges as 



lim^ 



^0+ r 



if IjI > 1. Thus, Bj must be zero for 



\j\ > 1 , and the condition for the occurrence of a singular 
solution is |j| < 1. So, the presence of an irregular solu- 
tion stems from the fact the operator is not self-adjoint 
for \j\ < 1, and this irregular solution is associated with 
a self-adjoint extension of the operator Hq [35, 36]. In 
other words, the self-adjoint extension essentially consists 
in including irregular solutions in V{Ho), which allows us 
to select an appropriate boundary condition for the prob- 
lem. 

In the present system the energy of a bound state has 



to be negative, so that k is a pure imaginary 
Thus, with the substitution k 



A, 



r\j\ 



iFi 



+Bj e-- p' 



^2 



in, we have {p' ~ 

b1,i + 2|j|,p') 

7-bU-2b-|,/ 



= IK. 

-2Kr) 



(17) 
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For Eq. (17) representing a bound state, the solution 
fij{r) must vanish for r — > oo, i.e., it must be nor- 
malizablc. By using the asymptotic representation of 
iFi{a,b, z) for r — > oo, the normahzabihty condition 
yields the relation 



i r(i + |j|) , 
r(i-b1) 



(18) 



From Eq. (16), for |j| < 1 we have Bj = Aj(-2k)2|jIAj; 
and by using Eq. (18), the bound state energy is 



£ = -- 



M 



1 r(i + bl) 



1 i/iii 



A,r(i-bl) 



(19) 



This coincides with Eq. (3.13) of Ref. [22] for a = 1, 
i.e., the spin- 1/2 AB problem in Euclidean space with 
the spinorial connection. By comparing Eq. (19) with 
Eq. (13), we find 



1 

a" 



1 



1 



m 

2 



2|j| 



1 



± m 



(20) 



We have thus attained a relation between the self-adjoint 
extension parameter and the physical parameters of the 
problem, j and rg. It should be mentioned that some 
relations involving the self-adjoint extension parameter 
and the (5-function coupling constant were previously ob- 
tained by using Green's function in Ref. [23] and the 
renormalization technique in Ref. [21], being both, how- 
ever, deprived from a clear physical interpretation. 

Once the bound energy problem has been examined, let 
us now analyze the AB scattering scenario. In this case, 
the boundary condition is again given by Eq. (14) but 
now with the replacement \j — ?> Aj, where A| is the self- 
adjoint extension parameter for the scattering problem. 
In the scattering analysis it is more convenient to use 
the solution of the equation "Ho/ij = k'^fij, in terms of 
Bessel functions 



hAr)=C,J{\j\,kr)+D,Y{\j\,kr), 



(21) 



with Cj and Dj being constants. Upon replacing fij{r) 
in the boundary condition (14), we obtain A^Cj^fc'-'l = 



Dj[Ck-\^\ 



A|(77fclj| C7fc"'^''linir->o+^^"^'^'')], where 



2ijir(i+b-|) 



. _ 2i^ir(|j|) _ 



cos(7r|j|)r(-|j|) 
7r2lJ 



and 
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,^■1 ' 1 .1, . As in the bound state calculation, when- 
ever \j\ < 1, we have Dj ^ 0; this means that there 
arises again the contribution of the irregular solution Y 
at the origin when the operator is not self-adjoint. Thus, 
for \2\ < 1, we obtain AJQCfcl^l = D^iCk'^^^ - AjT^fcl^l), 
and by substituting the values of ^, C a-nd rj into above 

expression we find Dj = —fi.'Cj, where 



A-fc^l^lr(l-ljl)sin(lj» 



A|pi^ir(i - \3\) cos(^IjI) + 4b-ir(i + \]\) ■ 



(22) 



Since the 5 is a short range potential, it follows that the 
behavior of /i j- for r — > oo is given by 



Trkr 



■ cos kr ^ 



TT 
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(23) 



where Sj^{k,(j)) is a scattering phase shift. The phase 
shift is a measure of the argument difference to the 
asymptotic behavior of the solution J(]m],fcr) of the 
radial- free equation that is regular at the origin. By us- 
ing the asymptotic behavior of J{\j\,kr) and y(]j],fcr) 
for r — > oo in Eq. (21), and comparing it with Eq. 

(23), similarly as done in [37], wc found that 5j^ (fc, cj)) ~ 

A„t{(j>) + where A„((/)) = f dm] - \m + (j>\), and 

9x^. = arctan(^^^). Therefore, the expression for the S 
matrix is 



,2iA„(0) 



A5pi,ir(i_|^-|)g,b>^4i,ir(i + |j|) 



A^fc2i^ir(i - ljl)e-'ij> + 4ijir(i + 



(24) 

In accordance with the general theory of scattering, the 
poles of the S matrix in the upper half of the complex 
plane [38] [these poles occur in the denominator of (24) 
with the replacement k — > in] determines the positions of 
the bound states in the energy scale, Eq. (19). From this, 
we have A| = Aj, with A^ given by Eq. (20), and the self- 
adjoint extension parameter for the scattering scenario 
being the same as that for the bound state problem. This 
is a very interesting result that has not been described 
in the literature yet, as far as we know. Thus, we also 
obtain the phase shift and the scattering matrix in terms 
of the physics of the problem. If A^ = 0, wc achieve the 
corresponding result for the pure AB problem with the 
Dirichlet boundary condition; in this case, we recover 
the expression for the scattering matrix found in Ref. 



„2»A„(0)^ If we make A? 



we get S = 



[39], S = e~ 

g2iA™(^)-|-2i7r|j 

In this article, we have presented a general regulariza- 
tion method to address a system endowed with a singular 
Hamiltonian (due to localized fields sources or quantum 
confinement). Using the KS approach, the bound states 
were determined in terms of the physics of the problem, 
in a very consistent way and without any arbitrary pa- 
rameter. In sequel, wc employed the BG approach; by 
comparing the results of these approaches, we have deter- 
mined the value of the self-adjoint extension parameter 
for the bound state problem, which coincides with the 
one for scattering problem. We thus obtain the S matrix 
in terms of the physics of the problem, as well. A natu- 
ral extension of the problem studied here, amongst many 
possible options, is the inclusion of the Coulomb poten- 
tial, which naturally appears in two-dimensional systems, 
such as graphene [40] and anyons systems [41, 42]. Re- 
sults in this respect will be reported elsewhere. 
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